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Abstract 

By exploiting the chiral symmetry, we derive analytic formulas for 
the massless and the massive Ginsparg- Wilson fermion propagators. Us- 
ing these formulas, we derive an expression for the chiral condensate 
which is the order parameter for spontaneous chiral symmetry break- 
ing in QCD. These formulas provide the proper way to compute the 
fermion propagators and may save a significant amount of computing 
time especially for large lattices in four dimensions. 
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It is well understood that the proper way to break the chiral symmetry of 
a chirally symmetric Dirac operator D c on the lattice is to replace D c by 

D = D C (1 + RD C )~ 1 (1) 

where R is any hermitian operator which is local in the position space and 
trivial in the Dirac space. Then the Dirac operator D satisfies the Ginsparg- 
Wilson relation |1| 

^75 + 75^ = 2D l5 RD (2) 

where the particular form of the chiral symmetry breaking on the RHS having 
one 75 sandwiched by two Dirac operators is the signature of the Ginsparg- 
Wilson fermion. This implies the chiral Ward identities, non-renormalization 
of vector and flavor non-singlet axial vector currents, and non-mixing of op- 
erators in different chiral representations, as shown by Hasenfratz ||. The 
general solution [Q, [3J] of the GW relation, Eq. ([!]), can be regarded as a chiral 
symmetry breaking transformation which gaurantees that a zero mode of D c 
is also a zero mode of D and vice versa, hence the index of D is equal to the 
index of D c . Then it is obvious to see that the presence of R cannot produce 
a zero mode for D if D c does not possess any zero modes in topologically non- 
trivial background gauge fields. We refer to ref. j|, |5| for further discussions 
on topological characteristics of D ( D c ). Now we can replace R in Eq. (0) 
by any two hermitian operators S and T having the same properties of R, and 
rewrite Eq. (|2]) as the following, 

£75(1 - SD) + (1 - DT) l5 D = , (3) 

then it is evident that the action A = -tpDip is invariant under the finite chiral 
transformation on the lattice 

^->^exp[«0(I-Dr)7 B ] (4) 
^exp[z07 5 (l-££>)]^ (5) 

where 9 is a global parameter. This generalizes the infinitesimal chiral transfor- 
mation observed by Luscher . In general, we can regard the Ginsparg- Wilson 
relation in the form of Eq. (|3]) as the exact chiral symmetry on the lattice, 
which recovers the usual chiral symmetry 75D + D75 = in the continuum 
limit ( a — > ). Such a generalization seems to be cruical for formulating 
chiral gauge theories on the lattice. 

In this paper, we will concentrate on the problem of evaluating the fermion 
propagator D~ l of the Ginsparg- Wilson fermion. From Eq. (^), we immedi- 
ately obtain 

IT 1 = LT 1 + R (6) 
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Hence, the fermion propagator D^ 1 is completely determined for any R if the 
chiral fermion propagator D^ 1 has been evaluated. In fact, the long distance 
properties of D^ 1 are governed by the long distance properties of D~ x since 
R is supposed to be local. Furthermore, we will show that the short dis- 
tance behavior of D -1 which is relevant to the the chiral condensate is also 
determined by D^ 1 . Therefore the physics of any Ginsparg- Wilson fermion 
operator D actually lies in its chiral limit D c . However, only D c itself is not 
sufficent to constitute the correct lattice fermion operator. According to the 
Nielson-Ninomiya no-go theorem @, if we require D to be local and free of 
species doubling, then the chiral symmetry {75, D} = must be broken. The 
transformation in Eq. ([!]) provides the proper chiral symmetry breaking which 
not only preserves the correct physics of D c , but also removes its non-locality 
for any gauge configurations as well as its singularities in topologically non- 
trivial gauge fields. We note in passing that the locality requirement for the 
GW fermion should be imposed on the unitary operator V in Eq. (|8|) rather 
than on D, since the non-locality of V is due to the initial construction and 
thus cannot be diminished, while the locality of D is always gauranteed by a 
properly chosen R in Eq. (JJ). 

In the following, we derive a formula for the chiral fermion propagator for 
any D c satisfying the hermiticity condition 



75J D c75 = D\ 



(7) 



The hermiticity condition and the chiral symmetry of D c implies that D c is 
antihermitian, thus there exists one to one correspondence between D c and a 
unitary operator V, 



D c = (l + V)(l-V)-\ V = (D C -1)(D C + 1) 



-1 



(8) 



where V also satisfies the hermiticity condition ( 75^75 = V* ). Then the 
unitary operator V can be expressed in terms of a hermitian operator h, 



V = l5 h 




hi hi 



hi 



h\ h 3 I -hi 




(9) 



where h\ = hi and h\ = /13. Using the unitarity condition V^V = J, we have 
h 2 = 1, 



h 2 



h\ + h 2 h\ h x h 2 + h 2 h 3 
hXhi + /13/4 h\h 2 + hi 



Then we immediately obtain 
D c = (l + V)(l-V)- 1 





-h\{l-hi 




[l-h^h 



,-1 



(10) 



11: 
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and 



d: 1 = {i-v){i + v)- 



o -(i + h)- 1 ^ 

^(1 + Zii)- 1 



(12) 



The last formula together with Eq. (||) reduce the task of computing D~ x for 
any R to that of computing (1 + h{)~ 1 only. For a four dimensional lattice of 
volume Q = L 4 , the former is to compute the inverse of a complex matrix of size 
4L 4 x 4L 4 ( multiplicities due to the color and the flavor have been suppressed 
), while the latter only amounts to computing the inverse of a hermitian matrix 
of size 2L 4 x 2L 4 . Further simplifications are possible, however, depend on the 
actual form of V ( D c ) as well as the algorithm used for computing the inverse. 
We do not intend to go into these technical details in this paper. In principle, 
the formula (|f~2|) provides the proper way to compute the fermion propagator 
for any Ginsparg- Wilson fermion and serves as a starting point for devising 
highly optimized algorithms for computing the fermion propagators. 

In order to investigate the spontaneous chiral symmetry breaking ( x$B ) 
of QCD, we usually introduce a mass parameter m to the Dirac operator D, 
and then measure the chiral condensate in the thermodynamic limit Q — > oo 
followed by the massless limit m — > 0. The proper way to introduce the 
mass parameter to the Dirac operator D is to add the mass m to the chirally 
symmetric D c as in the continuum theory, but not to D directly. Then the 
GW-Dirac operator (P becomes 

D = (D c + m)(l + RDcY 1 (13) 

It should be noted that the mass m is not added to the D c in the last inverse 
operator which is inserted for breaking the chiral symmetry a la Ginsparg- 
Wilson. Then the fermion propagator for massive GW fermion is 

IT 1 = R+ (l-mR)(D c + m)- 1 (14) 



where 



and 



n-i ( mCl-h^B -Bh 2 \ 

(Dc + mr= ( h\B m-Ml-/4fl(l-/n)-%] J (15) 



B = [l + m 2 + (l-m 2 )^}- 1 . (16) 



Again, the amount of computations using above formulas is significantly less 
than that of computing (D c + m)~ l directly. Once (D c + m)~ l is evaluated, 
then it almost takes no time to obtain D~ x for any R using Eq. ([14]). 
Now we can write down the usual expression for the chiral condensate 

X = lim lim - ( ^ tr \d~\x, x) 1) (17) 

m->0 Si^oo \l \ L R=0i / 
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where the brackets < • • ■ > denote averaging over each gauge configuration 
with the weight det(D) exp(— A g ), and A g is the pure gauge action. The reason 
to set R = in D~ x is obvious, since it is an arbitrary operator introduced 
for breaking the chiral symmetry a la Ginsparg- Wilson, thus it is irrelevant to 
the physical xSB of QCD, and must be dropped Q Inserting ( fUj ) into (O), 
we have 

X = lim lim U^2ti[(D c + m)-\x,x)}) (18) 



The subtleness of Eq. ( JJ~S| ) lies in its two limiting processes. If the order of these 
two limiting processes is exchanged, then \ must be zero since spontaneous 
chiral symmetry breaking could not occur on a finite lattice for exactly massless 
fermion. Using (|15l), we can simplify Eq. (|T8| ) to the following expression 

X = lim lim - (Y 2m tr [(1 - h x )B] (x, x) ) (19) 
m^o n^oo \ l y-^ I 

where B is defined in Eq. (pl|). This is a simpler expression than (|l^) since 
it only involves hi, the hermitian submatrix of h = 75 y defined in (^). The 
most time consuming operation in (|l9|) is to evaluate the inverse matrix B. 

Now we restore the weight det(D) exp(— A g ) embedded inside the brackets 
< • • • > and rewrite ([19]) in the following, 

1 1 

X 



where 



lim lim -- f[dU] exp(-A fl ) det(,D)V 2m tr[(l - h^BUx, x) (20) 

Z = J [dU] exp(-A g ) det(D) (21) 
and 

det(D) = det(D c ) det(l + RD c )~ l (22) 

using Eq. (|l]). From the GW relation (Q), it is obvious that we can redefine 
D' = RD such that D' satisfies (0) with R = 1. However, one cannot set R = 1 
to the valence quarks since that would produce terms which are not relevent 
to the physical chiral symmetry breaking in QCD, as discussed after Eq. (|17|). 
On the other hand, for quarks in the internal quark loops, the chiral symmetry 
must be broken such that D can be local and free of species doubling, then 
the fermionic determinant and the axial anomaly can turn out to be correct. 
Therefore we can set R to one in Eq. (^2|), provided that D c has the correct 
classical continuum limit. Using Eq. (|TTD, we obtain 

det(D) = det(i±^Udet(i^) (23) 



1 Note that setting R = in ([H]) does not necessarily reinstall doubling if D c is non-local, 
as have been discussed in detail in ref. j|. 
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In topologically nontrivial gauge sectors, the zero modes of D corresponding 
to the —1 eigenmodes of h\, but the +1 eigenmodes of h^. Further implemen- 
tations of computing x depend on the actual form of V ( or D c ) as well as the 
simulation algorithms. However, we will not go into these technical details in 
this paper. Now we put the results together and rewrite the chiral condensate 
in the following formula 



X 



1 + hx 



lim lim —\- [ \dU] exp(— A a ) det f 

m->0 fi-»oo ft Z ./ \ 2 J 

J2 2mti\(l-h 1 )[l + m 2 + (l-m 2 )h 1 }' 1 (x,x)] (24) 



where hi is the hermitian submatrix defined in Eq. (^]) and 

'1 + hi 



Z = J [dU] exp(-A 9 ) det (25) 



If the chirally symmetric D c is initially specified as 

-D\ ' 



D >={k j < 26 > 

but not in terms of V = 75/1, then it is straightforward to rewrite all formulas 
for the fermion propagators and the chiral condensate in terms of D^. 

Finally, we note that Neuberger |8| recently obtained an expression for the 
chiral condensate in the massless limit and in a fixed gauge background, in the 
framework of the overlap formalism ||. Neuberger's result would agree with 
our general formula (|24| ) if the mass parameter is introduced to the fermion 
propagator and the formula (|15|) is used. 

In summary, we have derived analytic formulas for the GW fermion prop- 
agators, for the massless and the massive cases respectively. These formulas 
provide the proper way to compute GW fermion propagators in analytic stud- 
ies as well as numerical ones. They also serve as a starting point for devising 
optimized algorithms for computing fermion propagators. An analytic formula 
for the chiral condensate is derived and is ready for numerical implementations 
for Monte Carlo simulations including dynamical fermions. 
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